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ABSTRACT 

We study the density profile and shape of the Galactic halo using deep multicolour images from the MENeaCS and CCCP projects, 
over 33 fields selected to avoid overlap with the Galactic plane. Using multicolour selection and PSF homogenization techniques we 
obtain catalogues of F stars (near-main sequence turnoff stars) out to Galactocentric distances up to 60kpc. Grouping nearby lines 
of sight, we construct the stellar density profiles through the halo in eight different directions by means of photometric parallaxes. 
Smooth halo models are then fitted to these profiles. We find clear evidence for a steepening of the density profile power law index 
around R = 20 kpc, from -2.50 ± 0.04 to -4.85 ± 0.04, and for a flattening of the halo towards the poles with best-fit axis ratio 
0.63 ± 0.02. Furthermore, we cannot rule out a mild triaxiality (w > 0.8). We recover the signatures of well-known substructure and 
streams that intersect our lines of sight. These results are consistent with those derived from wider but shallower surveys, and augur 
well for upcoming, wide-held surveys of comparable depth to our pencil beam surveys. 
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1. Introduction 


The stellar halo of the Milky Way only contains a tiny frac¬ 
tion of its stars, yet it provides important clues about the for¬ 
mation of the Galaxy and galaxy formation in general. Within 
the paradigm of hierarchical structure formation, galaxies evolve 
over time, growing by means of mergers and accretion of smaller 
systems. While in the central parts of galaxies the signatures of 
such events are rapidly dissipated, the long dynamical timescales 
allow accretion-induced substructures to linger for Gigayears in 
their outermost regions. Thus, the stellar structure of the outer 
halos of galaxies such as the Milky Way can help constrain not 
only the formation history of individual galaxies, but also cos¬ 
mological models of structure formation. 

Owing to the intrinsic faintness of stellar halos, the Milky 
Way is our best bet for a detailed study of such structures. How¬ 
ever, even studying the Galactic stellar halo is fraught with dif¬ 
ficulties; very sensitive data are required to probe stars at these 
large distances (out to 100 kpc), and spread over sufficiently 
large areas to constrain the overall structure as well as local¬ 
ized substructures. In recent decades the advent of CCD-based 
all-sky surveys such as the Slo an Digital Sky Survey (SDSS 
I York et al.l2000l:lAhn et al.l2014l) in the optical a nd the 2 Micron 
All Sky Survey (2MASS ISkrutskie et al.ll2006l) in the infrared 
have unlocked unprecedented views of the outer regions of the 
Galaxy. This has led to t he discovery of many previously un¬ 
known substructures (e . g.lNew berg et al. 2002; Belokurov et al. 


known substructures (e.g. ,\ewhere et al. ZUUZ; Belokurov et al 
2006b; Grillma^ l2006t (Belokuro v et al.ll2007l: lJuric et al.l l2008 


Bell et al.ll2008l) and to impr oved knowledge of the overall struc¬ 

ture in these outskirts (e.g . IChen et~aL||200 It Juric et ail 2008; 


Ide Jong etal.l20lot Sesar et al.l2010l.l201 lllFaccioli et al.l2014l) . 


Nevertheless, most of these recent analyses are still limited to 
either the inner parts of the stellar halo (Rcc ^ 30 kpc) or to 
particular, sparse stellar tracers (e.g. K-giants or RR Lyrae). 

In this paper we use deep photometry obtained with the 
Canada-France-Hawaii Telescope (CFHT) MegaCam and the 
Wide Field Camera (WFC) at the Isaac Newton Telescope (INT), 
scattered over a large range of Galactic latitudes and longitudes 
to probe main sequence turn-off (MSTO) stars out to distances 
of 60 kpc. Combining our data into eight independent lines of 
sight through the Galactic halo, we are able to constrain the over¬ 
all structure of the outer halo, and to probe the substructure in 
these outermost regions. In section 2 we describe the data set 
used for this analysis and the construction of our deep star cata¬ 
logues. Section 3 presents the derived stellar density profiles and 
smooth Galactic model fits. We discuss our results in section 4 
and present our conclusions in section 5. 


2. Observations and data processing 

2.1. Survey and observations 

We u se g and r images from the MENeaCS and the CCCP sur- 

veys (iSand et al .1120121 [Hoekstra et al.il201 2k iBildfelfe t al. 2012) 

together with several archival cluster fields from the CFHT- 
MegaCam instrument. We combine these data with U and i im¬ 
ages from a follow-up campaign with the INT-WFC instrument 
(van der Burg et al., in prep.). Whereas these surveys targeted a 
preselected sample of galaxy clusters, the pointings constitute a 
"blind" survey of the Milky Way stellar halo since their distri¬ 
bution is completely independent of any prior knowledge of the 
halo’s structure and substructure. 
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Fig. 1: Equatorial map showing the position of all the fields used in this work. The different colours and symbols indicate how the 
fields have been grouped to calculate the different density profiles. The background image is the SDSS-DR8 map from lKonosov et akl 
(120121) . which shows the footprint of the Sagittarius stream and the location of the Sagittarius dwarf galaxy. When grouping the 
fields, we have also taken into account the presence of this stream, the Triangulum-Andromeda overdensity, and the anticentre 
substructures (ACS, EBS, and Monoceros), in trying to combine their effect in certain profiles and avoid it in others. 


Our pointings are distributed over the region of the sky visi¬ 
ble to both the CFHT and the INT (see FigureQ}. To optimize the 
star-galaxy separation (see section [2~2l i we restrict our analysis to 
exposures with image quality of subarcsecond seeing, typically 
<~ 0.9 arcsec in the r band. This limitation, combined with the 
varying fields of view and observing conditions between the data 
sets, leads to pointing footprint sizes that range between 0.24 and 
1.14 deg 2 . 

2.2. Image correction of the PSF distortion [and implications 
for the star-galaxy separation] 

Previous research by our group has shown that the performance 
of standard star-galaxy separation methods based on the size and 
ellipticity of the sources can be improved by homogenizing the 
point-spread function (PSF)_across an image prior to its photo¬ 
metric analysis (iPila-Dfez et al.l 120141) . In addition, such a cor¬ 
rection also provides the benefit of allowing us to perform fixed 
aperture photometry and colour measurements. 

In order to homogenize the PSF of our images, we use a code 
JPila-Dfez et al.l[20l41) that, as a first step, takes the shape of the 
bright stars in a given image and uses it to map the varying PSF 
and, as a second step, convolves this map with a spatially vari¬ 
able kernel designed to transform everywhere the original PSF 
into a gaussian PSF. 

2.3. Catalogues 

From the PSF-homogenized expo sures we create photom etric 
catalogues using Source Extractor (iBertin & Arnoutslll996i) . For 
the g and the r data, we stack the different exposures in each 
band to create a single calibrated image, and we extract the 
band catalogues from them. We perform a star-galaxy separation 
based on the brightness, size and ellipticity of the sources and we 
match the surviving sources in the two catalogues to produce a 
gr-catalogue of stars for each field of view (see IPila-Dfez et ail 


(120141) '). The limiting magnitudes of these gr star catalogues 
reach m^B ~ 25.0 at the 5.0cr level in the r band. 

For the U and the i fields of view, we produce several photo¬ 
metric catalogues, one for each individual exposure. We correct 
the magnitudes in the i catalogues for the dependency of the il¬ 
lumination on pixel position. For each pointing and band, the 
exposure catalogues are calibrated to a common zero point and 
combined to produce a single-band catalogue. In these single¬ 
band catalogues, the resulting magnitude for each source is cal¬ 
culated as the median of the contributions of all the individual 
exposures. At this point the U and the i magnitudes are con¬ 
verted from the INT to the CFHT photometric system using the 
following equations, which we derive by calibrating our mixed 
INT-CFHT colours to the colour stellar loci of the CFHT Legacy 
Survey dErben et alJ d2009l) . [Hi Idehrandt et al.l d2009l) ): 

IMegaCam = IlNT ~ 0.12 * (rMega ~ ilNT) (1) 

U MegaCam = UlNT ~ 0.15 * (ujnT ~ gMega) ■ (2) 

Finally we position-match the sources from the U -, the i- and 
the g/--catalogues to create a final catalogue of stellar sources for 
each field of view. These final wgri-catalogues are shallower than 
the #/--catalogues because of the lesser depth of the i and the U 
observations ( see TableQ}. FigureQshows the colour-magnitude 
diagrams (CMDs) for the final ugri and gr catalogues (top and 
centre, respectively), and the difference between them (bottom). 
The bottom panel highlights that, in the colour regime of the halo 
(0.2 < g - r < 0.3), the combination of the four bands removes 
mainly very faint, unresolved galaxies. 

We correct for interstellar extinction using the maps from 
ISchlegel et al.l d 1998b and transform the magnitudes in the ugri- 
stellar catalogues from the CFHT to the SDSS photometric sys¬ 
tem. For this we use the equations on the Canadian Astronomy 
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Data Center MegaCam website 1 

U MegaCam = USDSS ~ 0.241 • (lisDSS ~ gSDSs) (3) 

g MegaCam = gSDSS ~ 0.153 • (gSDSS ~ r SDSS ) (4) 

fMegaCam = rgDSS ~ 0.024 • (gSDSS ~ ^SDSS ) (5) 

iMegaCam = h DSS ~ 0.003 ■ (r$ DSS ~ h DSs) (6) 


and invert them to turn our measurements into SDSS magni¬ 
tudes. Subsequently we calibrate each field directly to SDSS 
using stellar photometry from DR8. The resulting photometry 
matches the colour-colour stellar loci of ICovev et al. ( 2007 ) as 
shown in Figure[3] Unless explicitly stated otherwise, all magni¬ 
tudes in this paper are expressed in the SDSS system. 

In order to reduce the noise when analysing the radial stellar 
density distribution of the halo, we combine the catalogues from 
nearby pointings, grouping them according to their position in 
the sky. This step is important because of the nature of our sur¬ 
vey, which is composed of relatively small, scattered fields of 
view. We use a friends-of-friends (FoF) algorithm to group the 
different pointings. We request two friends not to be apart by 
more than 20 degrees, and in a few cases we clean or split a re¬ 
sulting group (red pentagons or blue and orange triangles in Fig¬ 
ure |T} or combine others (purple diamonds) to account for the 
positions of the galactic disk or major halo substructures. Be¬ 
cause the different pointings in our surveys have different com¬ 
pleteness limits, these grouped or combined catalogues -which 
we name A,B,C,... H- are finally filtered to meet the complete¬ 
ness magnitude threshold of their most restrictive contributor 2 . 


3. Stellar radial density profiles 

3. 1. Star selection and construction of the radial stellar 
density profiles 

The coordinates and the completeness limits of the groups are 
given in Table Q] We use halo main sequence turnoff stars in 
our fields as tracer of the stellar halo: at the completeness lim¬ 
its of the data such stars can be identified as far out as 60 kpc 
from the Galactic centre. We fit several Galactic stellar distri¬ 
bution models to these density profiles and derive a number of 
structural parameters for the stellar halo. Previous works have 
already used main sequence turnoff point (MSTO) stars, near- 
MSTO stars, BHB and blue stragglers of type A and RRLyrae 
as stellar tracers for the Galactic stellar halo. We compare and 
discuss our findings to theirs in section [4~2l 

In order to select the near main sequence turnoff stars we 
make use of two empirical photometric variables. The ratio 
[Fe/H] is calculated following the photometric metallicity re¬ 
lation by iBond et ahl ( 2010 ). and the absolute magnitude M r 
is calculated following the photometric parallax relation from 
llvezic et all d2008l) : 


[Fe/H] = -13.13+ 14.09x + 28.04y-5.51xy-5.90;r 

- 58.68y 2 + 9.14 Jy - 20.61xy 2 + 58.20y 3 , (7) 

1 www2.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/megapipe/docs/filters.html 

2 To determine the completeness limit of each field of view, we fit its 

magnitude distribution to a gaussian -representing the population of 
faint galaxies- and another variable function -representing the stellar 
distribution along the whole magnitude range-. We choose as the com¬ 
pleteness limit either the transition point between the two distributions 

(the valley) or, if instead there is a plateau, the turning point of the whole 

distribution (the knee). 


Distribution of sources from the griu catalogue 



Distribution of sources from the gr-only catalogue 
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Fig. 2: Hess diagrams showing the number of sources per colour- 
magnitude bin in the ugri catalogue (top), in the gr catalogue 
(centre) and the difference between both (bottom) for field 
A1033. Most of the sources lost when combining the catalogues 
correspond to faint magnitudes, because the i and the U obser¬ 
vations are shallower. The effect is the removal of most of the 
faint galaxies (located in the -0.2 < g - r < 0.7 and r > 23 
region in the central panel), most of the faintest disk M dwarves 
(1.1 < g - r < 1.3) and a number of faint objects (in the i or the 
U bands) scattered throughout the (g - r, r ) diagram. 


where x = u - g and y — g - r. This relation is valid in the 
g - i < 0.6 and -2.5 < [Fe/H] < 0 range, which is compatible 
with the regime of our near-MSTO star selection. 


M r = -0.56 + 14.32Z- 12.97z 2 + 6.127z 3 - 1.267z 4 

+ 0.0967z 5 - \.\\[Fe/H]-f).l%[FelH] 2 , (8) 
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Table 1: Groups of pointings as shown in Figures |T} [5] [6] and [8] The table shows the central coordinates for each group, the number 
of individual fields of view contributing to it, its total area and the stellar completeness limit in the r band. 


Group 

RA (deg) 

Dec (deg) 

l (deg) 

b (deg) 

ft fields 

2 (deg“) 

inagums,* (mag) 

A 

160.654338 

43.98310 

171.335811 

59.15040 

8 

5.60 

22.8 

B 

231.593130 

29.13513 

45.577138 

55.93598 

5 

3.98 

22.7 

C 

229.347757 

6.91624 

9.425402 

49.92775 

4 

3.44 

24.1 

D 

210.062933 

51.67173 

99.735627 

62.24580 

2 

0.64 

23.4 

E 

121.918411 

41.20348 

179.233500 

31.26694 

5 

2.73 

22.7 

F 

342.735895 

17.09581 

86.019738 

-36.99391 

3 

2.17 

23.2 

G 

157.028363 

17.15674 

222.142793 

55.48268 

3 

2.02 

23.1 

H 

220.659749 

2.00187 

354.337092 

53.38989 

3 

2.04 

24.2 



r - i 


Fig. 3: Colour-colour diagrams (CCDs) corresponding to the 
fields in group A (pointings marked as light cyan circles in Fig- 
ure|T|). The sources in the ugri catalogues (black) and the subset 
of near-MSTO stars (red) have been calibrated to SDSS using 
DR8 stellar photometry. The main sequence stellar loci (green 
dashed lines) are the ones given in Tables 3 and 4 of lCovev et all 
( 2007 1. Quasars and white dwarf-M dwarf pairs are abundant in 
the u - g < 1, -0.3 < g - r < 0.7 space. 


where z = g — i. The tested validity regime of this equation en¬ 
compasses the 0.2 < g — i < 1.0 range, meaning that the absolute 
brightnesses of our near-MSTO stars have been properly esti¬ 
mated. We extrapolate the relation for the 0.1 < g - i < 0.2 
range, which is justified owing to the smooth and slow change 
of M r with z. 



[Fe/H] 


Fig. 4: Estimated absolute magnitude in the r band (M r ) and esti¬ 
mated metallicity ([Fe/H]) for group A for the sources typically 
considered as halo stars (blue) and those that we have selected 
as near-MSTO stars (red). The sources selected as halo members 
meet 0.2 < g - r < 0.3 and g, r, i > 17. The subset of near- 
MSTO stars, additionally meets M r > -2, -2.5 < [Fe/H\ < 0 
and 0.1 < g - i < 0.6. 


We select the halo near-MSTO stars by requiring 


0.2 < g - r < 0.3 ; 

(9) 

g, r, i > 17 ; 

(10) 

0.1 < g - i < 0.6; 

(11) 

5.0 > M r > -2; 

(12) 

-2.5 < [Fe/H] < 0 . 

(13) 


The first two restrictions (l9landlT0l) retrieve stars typically as¬ 
sociated with the halo, in particular distant main sequence F stars 
(see Table 3 from lCovev et al. i ( 2007 )'). This selection however, 
can be significantly contaminated by quasars and white dwarf- 
M dwarf pairs, which are abundant in (but not restricted to) the 
-0.2 < g - r < 0.3 range (see Figure^. To reduce the presence 
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of these interlopers and select the bulk of the F stars population, 
we apply restrictions|TT](based on Table 4 in iCovev et ali d2007t) ) 
and[12] Constraint|T3]ensures that the final sources are at most as 
metal rich as the Sun (to account for possible contributions from 
metal-rich satellites) and not more metal-poor than 0.003 times 
the Sun. 

The decrease in interlopers attained by applying restric¬ 
tions [TT| HU and HU compared to only applying restrictions HI 
and HU is illustrated in Figure [3j where the red dots indicate the 
final selection of halo near-MSTO stars and the black dots rep¬ 
resent the whole catalogue of star-like sources. It is clear that the 
final selection of near-MSTO stars does not span the whole range 
of sources encompassed between g - r = 0.2 and g — r — 0.3. 
The effect of the [Fe/H] and M r selection is further illustrated 
in Figure |4] 

Using the estimated absolute brightness, we calculate the 
distance modulus and the heliocentric distance for all the near- 
MSTO stars. We define distance modulus bins of size Ap = 0.2 
mag and A// = 0.4 mag, and count the number of near-MSTO 
stars per bin for each group of fields (A,B,C,...). The choice of 
distance bins is motivated by a compromise between maximising 
the radial distance resolution and minimising the Poisson noise 
in the stellar number counts. We test this compromise by explor¬ 
ing two distance modulus bin sizes, which correspond to dis¬ 
tance bin sizes of the order of 10 2 pc and 10 kpc, respectively. 

We then calculate the number density per bin and its uncer¬ 
tainty as follows: 


Pl,b,D 

E p 


N IAAm 


0.2 ■ ln(l0) ■ D\ c ■ AQ-Ap 


Va" ~ fields 


(14) 

(15) 


Figure HI shows that the density profiles decrease quite 
smoothly for 40 - 60 kiloparsecs and for most of the lines of 
sight. 


3.2. Fitting procedure 

We fit several models of the Galactic stellar number density dis¬ 
tribution to the data, ranging from a basic axisymmetric power 
law to more complex models with triaxiality and a break in the 
power law. The models take the following mathematical forms, 
with x, y. and z being the cartesian galactocentric coordinates 
with the Sun at (8,0,0) kpc (lMalkinll2012l) : 


- Axisymmetric model 

p(x, y, z) = po • ^x 2 + y 2 + ^ j " /2 , (17) 

where q = c/a is the polar axis ratio or the oblateness of the 
halo; 

- Triaxial model 

p(x,y,z) = po • (x 2 + — + —\ " /2 , (18) 

\ w q) 


where w = b/a is the ratio between the axes in the Galactic 
plane; 

- Broken power law, with varying power index at Rbreak 


p(x,y,z) = 


f PO ■ 

(' R-ellip\ 

\ Po ■ 

( Rellip) 


n ou t . U n in~ n ou 
break 


Rellip ^ Rbreak 
Rellip — Rbreak 

z 2 


(19) 


Reilip — (x + y~ + 


1/2 . 


where AQ. is the area covered by each group. Dpi: is the helio¬ 
centric distance, / and b are the galactic coordinates and Nij,j\p 
is the number of stars per bin in a given direction of the sky. 
Particularly, 

4 JT 

AO = -X(deg 2 ) (16) 

41253 S 

and the area of each group (X) depends on the individual area of 
each field contributing to it (Table H}. 

The results for these number density calculations can be 
seen in Figure HJ where we plot the logarithmic number density 
against the galactocentric distance 3 , R^c , for each group (or line 
of sight). For this and the subsequent analysis, we only consider 
bins with Rcc > 5kpc, |z| > 10 kpc (to avoid the inner regions 
of the Galaxy) and a distance modulus of p < magu m - 4.5 (to 
guarantee a complete sample of the faintest near-MSTO stars 4 ). 


3 

Rgc = yjRac + z 2 

where Rcc and z are the radial and vertical coordinates on the cylindrical 
galactocentric reference system. 

4 This constraint guarantees that there are no distance completeness 
issues due to our specific type of stellar tracers and due to the dif¬ 
ferent depths of our fields. The only subset affected by incomplete¬ 
ness is that of magu m - 5.0 < p < magu m - 4.5 for the stars in the 
4.5 < M, < 5.0 range; and its average loss is of 20% over the total num¬ 
ber of near-MSTO stars (-2.0 < M r < 5.0) in the same distance range. 

Several tests on different upper distance thresholds for the density pro¬ 
files show that the distance modulus constraint of p < magu m - 4.5 is 


- Broken power law, with varying power index and oblateness 

at Rbreak 


p(x, y, z) 


PO ■ (x 2 + T + £) ”"‘ /2 , Rcc < Rbreak 
PO ■ (x 2 + y 2 + , Rcc > Rbreak , 

( 20 ) 


where the inner power law is fit to data that meets Rcc < 
Rbreak and the outer power law is applied to data that meets 
RCC > Rbreak- 


We fit all these models to the data using the "curve- 
fit" method from Scipy.optimize, which uses the Levenberg- 
Marquardt algorithm for non-linear least squares fitting. The ob¬ 
jective function takes the form of a % 2 , and we also calculate a 
reduced^- 2 for analysis purposes. 


A" 


2 

Xred 


Ndatci 

- I 

1=1 


Pdata,i Pmodel,i 


X 


Ndata - Ftparams 


( 21 ) 

( 22 ) 


where Ndata and Np ara ms are the number of data points and the 
number of free parameters, respectively. 

The influence of the photometric uncertainties on the den¬ 
sity profiles and the best fit parameters is evaluated through a 

enough to guarantee that all the lines of sight contribute robust density 
measurements at the furthest distances and that the incompleteness in 
magi im - 5.0 < p < magu m - 4.5 for the 4.5 < M r < 5.0 near-MSTO 
stars has no statistically significant effect on the best fit parameters. 
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Fig. 5: Logarithmic stellar density profiles versus distance for the near Main Sequence turnoff point stars (near-MSTO) from the 
fields in groups A (green circles), B (cyan squares), C (blue downward triangles), D (yellow upward triangles), E (red pentagons), 
F (pink hexagons), G (purple diamonds) and H (orange leftward triangles). Their symbols match those in Figure!]] 
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set of Monte Carlo simulations that randomly modify the g,r,i,u 
magnitudes of each star within the limits of the photometric un¬ 
certainties. Through this method we find that the variation of the 
Monte Carlo best fit parameters aligns with the uncertainties of 
our best fit parameters (derived from the second derivative of the 
fits by the "curve-fit" method). The centre of these variations is 
within lcr of our direct findings. 

We fit all models to four data sets: with and without [known] 
substructures and binned in 0.2 and 0.4 magnitude cells. In this 
way we can check the robustness of our results to different bin¬ 
ning options and we are able to compare what would be the effect 
of substructure on our understanding of the smooth halo if we 
were to ignore it or unable to recognize it as such. Specifically, 
we cut the distance bins at Rcc < 25 kpc in group E to avoid 
contributions by the structures in the direction of the galactic 
anticentre (the Monoceros ring, the Anticentre Structure and the 
Eastern Band Structure), the distance bins within 15 < D/,c < 40 
kpc in group G to avoid contributions by the Sagittarius stream, 
and the distance bins within 20 kpc < Z)/,c < 60 kpc in group H 
to avoid contributions again by the Sagittarius stream. 

3.3. Results 

The best fit parameters for each model resulting from fitting 
these four data sets are summarized in Tables [2] to [5] Table [2] 
contains the results of fitting the A/j = 0.2 mag binned data ex¬ 
cluding regions with substructure, whereas Table —contains the 
results of fitting to all the 0.2 mag bins. Similarly Table-cov¬ 
ers the fits to A/j = 0.4 mag data without substructure bins, and 
Table — to all 0.4 mag bins. The reduced^ 2 and the initial pa¬ 
rameters have also been recorded in these tables. 

We compare the fitting results for the four different data sets 
recorded in Tables —to —and find that the fits for which the sub¬ 
structure has been masked significantly outperform those that 
have been allowed to fit all the available data. The difference on 
X 2 m/ for all these models and bin sizes is in every case at least 
a factor of 2.3 or larger. We find that allowing the models to fit 
data that contains substructure does not affect largely most of the 
structural parameters (polar axis ratios are compatible within the 
uncertainties and power law indices have close values) except 
that it decreases the disk axis ratio w by at least 15%, suggest¬ 
ing a strong departure from the axisymmetric model that is not 
implicit in the filtered data sets. Henceforth we will restrict the 
remaining discussion to the results derived from the cleanest data 
sets. 

Comparing the parameters resulting from the best fits to the 
masked 0.2 mag and 0.4 mag data, we find that the fits to 0.2 
mag binned data perform better for all the models (x 2 d ratio of 
two). Nonetheless, all the measurements for the different struc¬ 
tural parameters in the two data sets are compatible with each 
other within the uncertainties. The best fits for the four models 
and their residuals for our eight lines of sight are shown in Fig¬ 
ures [6a] and [6b] for the masked 0.2 mag binned data. It is clear 
that the differences between the fitted models along these sight 
lines are small. 

Our data are inconclusive regarding triaxiality, but are com¬ 
patible with either a moderately triaxial halo or with no triaxial¬ 
ity. For the 0.2 mag data set, the triaxial model fits slightly better 
than the axisymmetric model and returns w = 0.75+0.09. For the 
0.4 mag data set, however, the axisymmetric model fits slightly 
better and the triaxial model returns a disk axis ratio compati¬ 
ble with 1. In both data sets the other best-fitting parameters are 
practically identical for the two models. This indicates that the 
cost of the extra parameter is not supported by the 0.4 mag data. 


Thus, it is hard to derive a precise value for the disk axis ratio 
and to conclude if it is truly triaxial, but a weighted average of w 
and the general analysis show confidently that w > 0.8. 

We increase the complexity of the axisymmetric model by 
adding two degrees of freedom and considering a change in the 
power law index n at a specific break distance Rbreak (a broken 
power law). For this purpose, we use a grid of values to ex¬ 
plore all the parameters except the density scale factor f>o , which 
we left free to fit (see below for the grid characterization). This 
model decreases the x 2 red in both the 0.2 and the 0.4 mag binned 
cases, indicating that our data is better fit by a broken power 
law than by a simple axisymmetric model or a triaxial model. 
It turns the single power law index from n = -4.26 ± 0.06 into 
a less steep inner index n,„ = -2.50 ± 0.04 and a steeper outer 
index n out = -4.85 + 0.04 (measurements here are for weighted 
averages between the 0.2 and 0.4 mag data). It also increases the 
central value of the polar axis ratio q within the uncertainties, 
from a weighted q = 0.60 ± 0.04 to a weighted q = 0.63 + 0.02. 
Globally, the disk axis ratio seems to be the most stable param¬ 
eter throughout the different model fits to our data, returning a 
quite oblate halo. 

Finally we fix the break distance at the best fit value found 
by the broken power law model (Rbreak = 19 kpc and 20 kpc for 
the 0.2 and 0.4 mag binned data, respectively) and add another 
parameter to it, allowing not only n, but also q to change at the 
break distance. We find that the best fits to this model return 
such large error bars for the inner halo that, in practice, it yields 
unconstrained measurements: Apy < po, An,-,, is 12-18% of 
and A q, n is 40% of q m . 

We explore each model to investigate possible parameter de¬ 
generacies, tolerance ranges and potential local minima in our 
best fits. For this we fix all the parameters in the four models 
except the density scale factor po, and we run the fits across a 
grid of parameter values. In particular, the grids are built fol¬ 
lowing q,w e [0.1,2.0;b = 0.05], n e [-5.0- 1.0;<5 = 0.1], 

e [-4.0,-1.0; 6 = 0.1], n oll , 6 [—7.0, —3.0; <5 = 0.2] and 
Rbr e [15,50; 6 = 1], where 6 is the incremental step for each 
parameter. We find that there is a degeneracy between Rb r and 
n in for the simple broken power law model for both binnings 
(see Figure—. 

Finally our measurements for the density scale factor po (p 
at Rcc = 1 kpc) are the result of large extrapolations and merely 
serve as normalizations for our fits. For that reason we do not 
discuss these values in detail. 


4. Discussion 

4. 1. Robustness of the best fit structural parameters 

In order to determine how the data available to us influences the 
results from our best fits, we remove the different lines of sight 
one at a time and repeat the fits. In this way we can determine 
which are the most critical lines of sight and what is their effect 
on our results. 

We find that most of them have no significant influence on 
the best fit parameters of the different halo models. However, 
starting with the polar axis ratio we find that removing group A 
significantly increases its value (q ~ 0.7) and removing groups C 
or E significantly decreases it (q ~ 0.5) in both the axisymmet¬ 
ric and triaxial model in the two data sets. Regarding the power 
law index, again groups A or C have an influence, but group B 
as well. Removing groups A or B increases n to ~ -4.1 ±0.1, 
whereas removing C decreases it to n ~ -4.6. When considering 
a triaxial halo, we find that groups A, B or C increase the disk 
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Table 2: Best fit parameters for the four different Galactic stellar distribution models resulting from removing the data that is affected by known halo substructures (the Sagittarius 
stream and the anticentre substructures). For the fitting, the data has been binned in 0.2 mag distance modulus cells. 


Model 

xL 

po(pc- J )-10- J 

Rbreak (kpc) 

n 

Min 

M-out 

q 

Qin 

tfout 

W 

axisymmetric 

1.90 

14 + 6 

- 

-4.31 + 0.09 

- 

- 

0.62 ± 0.06 

- 

- 

- 

triaxial 

1.86 

14 + 6 

- 

-4.28 ± 0.09 

- 

- 

0.60 ± 0.06 

- 

- 

0.75 ± 0.09 

broken p.l.„ 

1.52 

0.071+0.003 

19.0 + 0.5 

- 

-2.40 ± 0.05 

-4.80 + 0.05 

0.60 ± 0.03 

- 

- 

- 

broken p.l.„, 9 

1.99,1.51 

1+3 

19 fixed 

- 

-3.3 ±0.6 

-4.9 ± 0.2 

- 

0.5 ± 0.2 

0.68 ± 0.07 

- 

initial parameters 

- 

0.001 

40.0 

-3.00 

-3.00 

-3.50 

0.70 

0.70 

0.8 

1.00 


Table 3: Same as in Table[2]but this time fitting all the available data (including those regions containing stellar counts from known substructures and detected overdensities). 


Model 

xied 

Po (pc J ) ■ 10 J 

R-break (kpc) 

n 

Min 

Mout 

q 

Qin 

tfout 

W 

axisymmetric 

4.71 

8 + 3 

- 

-4.15 + 0.08 

- 

- 

0.69 ± 0.06 

- 

- 

- 

triaxial 

4.59 

7 + 2 

- 

-4.07 ± 0.08 

- 

- 

0.68 ± 0.06 

- 

- 

0.59 ± 0.07 

broken p.l.„ 

4.24 

0.17 + 0.01 

21.0 ±0.5 

- 

-2.80 + 0.05 

-4.80 + 0.05 

0.70 ± 0.03 

- 

- 

- 

broken p.l.„, 9 

3.36,4.79 

1+2 

21 fixed 

- 

-3.3 ±0.4 

-5.0 ±0.2 

- 

0.5 ± 0.2 

0.79 ± 0.08 

- 

initial parameters 

- 

0.001 

40.0 

-3.00 

-3.00 

-3.50 

0.70 

0.70 

0.8 

1.00 


Table 4: Same as in Table[2]but this time fitting the data binned in 0.4 mag distance modulus cells. 


Model 

Xfed 

po(pc^)-10- J 

Rbreak (kpc) 

n 

Min 

Mout 

q 

Qin 

tfout 

W 

axisymmetric 

3.89 

12 ±4 

- 

-4.26 ± 0.08 

- 

- 

0.60 ± 0.05 

- 

- 

- 

triaxial 

3.97 

12 ±5 

- 

-4.25 ± 0.08 

- 

- 

0.60 ± 0.06 

- 

- 

0.9 ±0.1 

broken p.l.„ 

2.61 

0.11+0.01 

20.0 ± 0.5 

- 

-2.60 + 0.05 

-4.90 ± 0.05 

0.65 ± 0.03 

- 

- 

- 

broken p.l.„ >? 

4.95,2.34 

1 ± 1 

20 fixed 

- 

-3.2 ±0.4 

-5.0 ±0.3 

- 

0.5 ± 0.2 

0.67 ± 0.08 

- 

initial parameters 

- 

0.001 

40.0 

-3.00 

-3.00 

-3.50 

0.70 

0.70 

0.8 

1.00 


Table 5: Same as in Table[4]but this time fitting all the available data (including those regions containing stellar counts from known substructures and detected overdensities). 


Model 

XL 

Po (pc j ) • 10 3 

Rbreak (kpc) 

n 

Min 

Mout 

q 

1 Jin 

tfout 

W 

axisymmetric 

9.13 

7 + 2 

- 

-4.10 + 0.07 

- 

- 

0.66 ± 0.05 

- 

- 

- 

triaxial 

9.19 

7 + 2 

- 

-4.07 ± 0.07 

- 

- 

0.65 ± 0.06 

- 

- 

0.74 ± 0.09 

broken p.l.„ 

7.74 

0.058 ± 0.005 

20.0 ± 0.05 

- 

-2.40 ± 0.05 

-4.8 ± 0.05 

0.70 + 0.03 

- 

- 

- 

broken p.l ., uq 

6.05,9.2 

0.6 ± 0.9 

20 fixed 

- 

-3.1 ±0.4 

-4.9 ± 0.2 

- 

0.5 ± 0.2 

0.74 ± 0.07 

- 

initial parameters 

- 

0.001 

40.0 

-3.00 

-3.00 

-3.50 

0.70 

0.70 

0.8 

1.00 
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(a) Fitted density profiles for the 0.2 mag binned data. 
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(b) Data-to-model residuals for the 0.2 mag binned data. 
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Fig. 6: Left panels: density profiles in decimal logarithmic scale and the best fit models from Table [2] (fitted to masked 0.2 binned data). Right panels: Residuals between the data 
and the best fit models from panel[6a] The different lines represent the axisymmetric (black solid line), the triaxial (green dashed line), the broken power law with varying power 
index (red dotted line) and the broken power law with varying power index and oblateness (blue dashed-dotted- dotted line) models. The grey areas denote data that have been 
masked from the fitting to account for the presence of substructure. 


B. Pila-Dfez et ah: A skewer survey of the Galactic halo from deep CFHT and INT images 












































































A&A proofs: manuscript no. MWstellarhalo_structure_PilaDiez 



(a)/tr e rf map for the filtered 0.2 mag binned data set. (b )xl d ma P f° r the filtered 0.4 mag binned data set. 


Fig. 7 '■ X 2 re d isocontours maps for and R/,,- from the simple broken power law model. The minimum is indicated with a white star. 
The black solid isocontours range from min(x 2 red ) + 0.1 to the maximum value, whereas the white dashed isocontours range from 
min(x 2 r .) + 0.01 to min(x 2 red ) + 0.05. The maps illustrate a degeneracy between both parameters in the best fits. 


axis ratio w by 0.10 - 0.15, and that removing groups E or F 
decreases it to w ~ 0.5. Additionally, in conditions of triaxiality, 
the lack of group E reduces q further to q » 0.35. 

Thus removing group E turns out to be critical for both q 
and w, representing a totally differently looking halo (extremely 
oblate and quite elliptical in the plane). Group F also has a sim¬ 
ilar effect on w but not on q. The reason why group E has such 
a strong influence in the determination of a possible triaxiality 
is that it is by far the closest group to the Galactic anticentre. 
Other groups also influence the measurements of the different 
parameters, but have a smaller influence on the general picture 
we would derive. Overall we see that the lines of sight we use 
can have a drastic effect on the w results and a significant but 
moderate effect on q and n. This means that a global view of the 
halo is essential owing to its complex structure. 

4.2. Comparison to previous studies 

Previous investigations using near-MSTO stars have explored 
both the inner and the outer halo out to moderate distances 
(30 - 40 kpc), and similar regimes have been probed with blue 
horizontal branch stars and blue struggler stars, MSTO stars or 
multiple stellar halo tracers. Studies involving RRLyrae stars 
have reached further out to 50kpc. Remarkably, the depth of our 
data allows us to probe further than any previous study (out to 
60 kpc) in several directions, independently of the stellar tracer. 

In this section we compare our findings regarding the struc¬ 
tural parameters of the stellar halo to those of the following re¬ 
sults in the literature: 

- iJuric et all d2008l) use near-MSTO stars from the SDSS-DR3 
and DR4 as stellar tracers, and cover the 5 kpc < Rcc <15 
kpc range. They comprise 5450 deg 2 in the northern Galactic 
hemisphere and 1088 deg 2 in the south. 

- ISesar et ahl ( 2011 ) use as well near-MSTO stars from the 
CFHT Legacy Survey, and explore the 5 kpc < Rcc < 35 
kpc range. Two of their four fields explore the South Galac¬ 
tic Cap. 


- iDeason et al.l (i201 ll) use type A blue horizontal branch 
(BHB) stars and blue stragglers (BS), reaching out to Rcc = 
40kpc. 

- Ide Jong et al.l (120101) use CMD fitting of SEGUE stellar pho¬ 
tometry to probe the total stellar mass density from Rcc = 7 
kpc to Rqc = 30 kpc along a "picket fence" of 2.5 degree 
wide strips at fixed Galactic longitude spanning a large range 
of Galactic latitudes. 

- IChen et akl (1200 ll) use more general MSTO stars from two 
high latitude regions of SDSS to the North and the South of 
the Galactic plane (49 deg < \b\ < 64 deg). They explore the 
inner halo re gime (Rgc ^ 30 kpc). 

- iBell et all (120081) use also more general MSTO stars from 
SDSS-DR5 spanning 5 < Rcc < 40 kpc. 

- iFaccioli et al . (1201 4) use RRLyrae in the 9 kpc < Rcc < 
49 kpc range. Their multiepoch data comes from the Xuyi 
Schmidt Telescope Photometric Survey (XSTPS) in combi¬ 
nation with SDSS colours, and covers 376.75 deg 2 at RA ~ 
150 deg and Dec « 27 deg. 

- ISesar et akl (1201 Ol) use RRLyrae stars from SDSS-II in the 
stripe 82 region. Although their data originally spans 5 k pc < 
Rcc <110 kpc, the reanalysis performed bv IFaccioli et al.l 
(20141) to derive structural parameters truncates the sample 
at 49 kpc. 

- IWatkins et al.l (l2009t) use as well RRLyrae from SDSS in 
stripe 82, and the comparative derivation of structural pa¬ 
rameters bv lFaccioli et al.l ( 2014 ) also truncates it at 49 kpc. 
Stripe 82 is located in the South Galactic Cap. 


The result of this comparison is summarized in Table [6] 
We note that the obla teness values for Faccioli et a l. (2014), 
ISesar et akl (1201 Ot) and IWatkins et al.1 (2009) are not the result 
of absolute best fits to a set of free parameters, but the best fits 
to free Rb r , n m and n out with fixed prior values for a quite oblate 
(q = 0.59^qq 2 ) and a moderately oblate halo (q = 0.70 ± 0.01). 

All surveys that reach beyond Rcc = 30 kpc coincide in 
the need for a break in the power-law index of the halo density. 
Regarding possible triaxiality, only a few of the studies report 
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constraints on w. Those that_do, have either reported ’finding 
unreasonable values’ (ISesar et al.ll201 1 1) or have obtained limits 
on triaxiality similar to ours (w > 0.8, Bell et al.ld2008h ’>. 

On the break radius, there is a general consensus towards 
Rbreak ~ 27 kpc. The only exception is that of iBell et alJ (i2008l) . 
who find a value very close to our measurement (~ 20 kpc). 
These discrepancies, however, can be explained by the effect of 
the Rbreak-nin degeneracy discussed in section [331 

The inner and outer halo power law indices mostly fall in 
the [-2.3,-3.0] and [-3.6,-5.1] ranges. Our inner power law 
index n,„ = -2.50 ± 0.04 is consistent with these results, par¬ 
ticularly with the lower end. In the case of the outer halo power 
inde x (n nu , = -4 .85 ± 0.04), the comparison is le ss trivial. First, 
only ISesar et ail (1201 ll) and iDeason et al.l (1201 lh have provided 
measurements for n out based on fits with a free q parameter 
{n out = -3.8 ± 0.1 and —4.6*[}p respectively). Second, only 
one work with n out measurements dSesar et al.ll201 lb uses a stel¬ 
lar tracer similar to ours (the others use A-BHB and BS stars, 
or RRLyrae stars). Most important, a good constraint on n out 
requires deep data, and none of these earlier surveys reach as 
deep as our data set. Our steep outer index, although well in the 
range of previous measurements, might well indicate a progres¬ 
sive steepening of the halo density, though it would be good to 
test this with additional sight lines of comparable depth. In any 
case, it seems safe to conclude that n out < -4.0. 

The best fit values for the polar axis ratio or oblateness q 
range from 0.5 to 0.9, with most of the measurements concen¬ 
trated within (0.55,0.70). The values of q do not seem to de¬ 
pend on whether a break was detected or not, nor on the limiting 
distance of the survey or on the stellar tracer. The discrepan¬ 
cies can thus be attributed either to methodological differences 
or to differences in the spatial coverage of the data samples. 
However, it is difficult to determine the actual cause. Our re¬ 
sults fit well within the most constricted range, and the simple 
broken power law measure (which is the one with the best x 2 r d 
and q = 0.63 + 0.04) fall s near t he upper end, bein g in good 
a greement with lJuric et all ([2008 h Sesar et akl (1201 lh and close 
to lDeason et al.l d20 1 lh and lChen et al.l ( 20011) . 

Finally it is noteworthy that the choice of stellar tracer across 
the different works does not seem to cause any significant bias 
on the best fit parameters. 


and in models dLaw & Maiewskil (l2010h and iPefiarrubia et al.l 
d20l(il) ). we expect this stream to span the 20 < D^c < 60 kpc 
or 16 < Rqc < 55 kpc range at these coordinates. This expecta¬ 
tion is met all along: they steadily increase from Rcc ~ 15 kpc, 
depart from p/pm = 3 ± 1 at Rcc = 30 kpc, reach p/p^r = 6 + 2 
at Rcc = 40 kpc and peak at Rqc = 45 kpc with max(p/pM ) = 
(12,15) ± 2. However, they do not decrease near R ( ;c = 55 kpc 
but seem to stay stable with a significant p/pm > 7+2). This sug¬ 
gests a thicker branch than predicted by the models, but in agree ¬ 
ment with prev i ous RRLyrae measurements (llbata et al] (1200 lh . 
iTotten & Irwinl (119981) and Pohm-PalmeretaiT i 200 lh as sum¬ 
marized in Figure 17 of lMaiewski et al. ( 2003h l . 

Two more modest overdensities that do not appear in the lit¬ 
erature seem to be present in groups C and D. In group C, a weak 
but consistent overdensity spans a distance range of Rcc ~ 35 
kpc to Rqc ~ 60 kpc. In group D, a sharp bump extends over a 
few kiloparsecs around Rcc < 20 kpc. 

We have looked for other known overdensities that position- 
match our lines of sight (see Figure [8}, but found no indication 
of them in the residuals. The first one corresponds to the tida l 
tails of the NGC5466 globular cluster (Belokurov et al.ll2006a f. 
which overlap with one field in group A and another one in group 
B (A1361 centred at (RA, Dec) = (176.09,46.39) and A1927 at 
(RA, Dec ) = (217.92,25.67)). This is a very weak cold substruc¬ 
ture located at Rcc ~= 16 kpc and extending for 45 d eg with an 
average width of 1.4 deg (IGrillmair & Johnsonll2006h . As such, 
it is not surprising to find no signature in the density profiles. 

The second one is the ensemble of three known overdensities 
in th e direction of group E: the Anti Center Stream (Rcc =18+2 
kpc, iRocha-Pinto et alJ d2003h and Li_eta[J (1201 2h l. the Mono- 
ceros ring (Rgc ~ 18 knc. lLi etak ( 201 2h ) and the Eastern Band 
Structure (Rcc = 20 + 2 knc. lLi et al.l ( 2012h l. These substruc¬ 
tures are masked from our fits and residuals when we impose 
|z| > 10 kpc to avoid the influence of the thick disk, and there¬ 
fore, they cannot be detected. 

The Triangulum-Andromeda overdensity ( (Martin et al.! 
12007 )1 falls close to one of the fields in group F. Despite this 
proximity, the residuals show no evidence for an overdensity at 
the expected distance of Rcc ~ 30 kpc, indicating that the over¬ 
density does not extend further in this direction. 


4.3. Detection of overdensities and identification 

We analyse the data-to-models residuals for the different lines 
of sight in Figure [6b] in search for overdensities. We find that, 
in general, all the lines of sight present regions with data-to- 
models deviations of a maximum factor of two. Additionally, 
certain lines of sight -C,D, G, and H- present more significant 
deviations spanning from a few kiloparsecs to tens of kilopar¬ 
secs in distance. We discuss these overdensities in greater detail 
below, and we also discuss expected overdensities that show no 
signature in our data. 

The most prominent overdensities in the data-to-model resid¬ 
uals correspond to the northern wrap of the Sagittarius (Sgr) 
stream. This stream overlaps in projection with groups G and H 
(see Figure [8}. For group G, the residuals indicate overdensities 
in the distance range where we expect to find both the Sgr and 
the Orphan strea m (2 0 < D^c < 40 kpc or 25 < Dec 5 44 kpc, 
iPila-Dfez et al] (l2014l lf. The overdensities indeed peak between 
Rcc = 25 kpc and 45 kpc, reaching p/pM =7 + 2, and drop 
sharply afterwards. Group H probes the Sgr stream closer to the 
Galactic centre but also for larger distances than group G. Based 
both on extensive data (summarized in iPila-Dfez et al . ( 1201 4-11 


5. Conclusions 

In this paper we have used wide-field images from the CFHT and 
the INT telescopes in eight broad lines of sight spread across 
the sky to produce deep photometric catalogues of halo near 
main sequence turnoff (near-MSTO) stars. Our images have been 
corrected for PSF inhomogeneities, resulting in catalogues with 
fixed-aperture colour measurements and improved star-galaxy 
separation. Thanks to the depth and quality of our data, we reach 
stellar completeness limits ranging from 22.7 mag to 24.2 mag 
in the r band, which translate into a 60 kpc distance limit for 
near-MSTO stars. 

We calculate galactocentric distances for the stars based on 
the photometric paralla x method by Ivezic et all (120081) and the 
metallicity estimator by Bon d et all ( 201 Oh . Webin them by dis¬ 
tance modulus, and calculate the stellar number density distribu¬ 
tion along the eight different lines of sight. 

In selecting the halo near-MSTO stars, we have used ad¬ 
ditional constraints than the standard 0.2 < g - r < 0.3 and 
g, r , i > 17 cuts in order to obtain a cleaner sample. Particularly, 
by applying additional cuts based on g-i colour, absolute mag¬ 
nitude and metallicity, we get a sample of mainly F stars signif- 
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Table 6: Comparison between the best fit structural parameters found in this work (weighted averages for the parameters of the 0.2 and 0.4 mag data sets) and those reported by 
other groups i n previous work s. The di fferent works ha ve bee n labelled as follows : J08 (lJuric et al. 20081) , SI 1 (ISesar et al. 20jl ), D11 ( Deas on et al.ll20111) . dJIO tide Jong et alJ 
1201 Oh . ChOl dChen et al.ll200H) . B08 dBell et al.ll2008h . F14 dFaccioli et al.ll2014h . and S10 (ISesar et al.ll20Toh and W09 dWatkins et al.ll2009h as reanalysed in F14. The fitted 
models in F14, S10 and W09 have fixed oblateness and test two different values motivated by the previous findings in SI 1 and Dll. 


Work 

stellar tracer 

dist. range (kpc) 


R br (kpc) 

n 

Win 

W-OUt 

q 

w 

this work-axisym. 

near-MSTO 

[10,60] 

1.9 

- 

-4.28 ± 0.06 

- 

- 

0.61 + 0.04 

- 

this work-triax. 

near-MSTO 

[10,60] 

1.9 

- 

-4.26 ± 0.06 

- 

- 

0.60 ± 0.04 

0.81+0.07 

this work-broken 

near-MSTO 

[10,60] 

1.5 

19.5 + 0.4 

- 

-2.50 + 0.04 

-4.85 ± 0.04 

0.63 ± 0.02 

- 

J08 

near-MSTO 

[5,15] 

[2,3] 

- 

- 

-2.8 ± 0.3 

- 

0.65 + 0.15 

- 

Sll 

near-MSTO 

[5,35] 

3.9 

27.8 + 0.8 

- 

-2.62 ± 0.04 

-3.8 ±0.1 

0.70 ± 0.02 

excluded 

Dll 

A-BHB, -BS 

[-,40] 

- 

27.1 ± 1 

- 

-2.3 ±0.1 

a fA- 0.2 

• u 0 1 

0 sq+O-O 2 

0 03 

- 

dJIO 

multiple 

[7,30] 

[3.9,4.2] 

- 

-2.75 ± 0.07 

- 


0.88 ± 0.03 

- 

ChOl 

MSTO 

[-,30] 

- 

- 

-2.5 ±0.3 

- 

- 

0.55 ± 0.06 

- 

B08 

MSTO 

[5,40] 

2.2 

~ 20 

-3 + 1 

- 

- 

[0.5,0.8] 

> 0.8 

F14 

RRLyrae 

[9,49] 

0.8 

28.5 + 5.6 

- 

-2.8 ± 0.4 

-4.4 ± 0.7 

q fix = 0.70 ± 0.01 

- 

tt 

RRLyrae 

[9,49] 

1.04 

26.5 + 8.9 

- 

-2.7 ± 0.6 

-3.6 ±0.4 

af - o 59 +0 - 02 
qfix u - J ^-o.o3 

- 

S10 

RRLyrae 

[9,49] 

1.1 

34.6 + 2.8 

- 

-2.8 ± 0.2 

-5.8 ±0.9 

q/ix = 0.70 ± o.oi 

- 

It 

RRLyrae 

[9,49] 

1.52 

26.2 ± 7.4 

- 

-3.0 ±0.3 

-3.8 ±0.3 

?/« = 0.59!» : “ 

- 

W09 

RRLyrae 

[9,49] 

1.1 

27.6 + 3.3 

- 

-2.5 ± 0.3 

-4.3 ± 0.4 

qfix = 0.70 ± 0.01 

- 

M 

RRLyrae 

[9,49] 

0.69 

26.9 + 3.1 

- 

-2.1 ±0.3 

-4.0 ±0.3 

q f i X = 0.59*% 

- 
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Fig. 8: Equatorial map showing the position of all the fields used in this work and the closest cold stellar overdensities to them. These 
overdensities are used for comparison and discussion of the stellar density profile data-to-model residuals throughout section l4~3l 
The labels in the figure correspond to the Anticentre Structure (ACS), the Eastern Band Structure (EBS), the NGC5466 stream, 
the Grillmair & Dionatos stream (G&D), the Orphan stream, the Triangulum-Andromeda overdensity (Tri-And) and the Pisces 
overdensity. The background image is the SDSS-DR8 map from lKonosov et ah ( 2012 ). which shows the footprint of the Sagittarius 
stream. The Monoceros ring also appears partially in this background image, as a dark region overlapping the western part of the 
Galactic disk in the anticentre region, eastwards of the ACS. 


icantly decontaminated from quasars and white dwarf-M dwarf 
pairs. 

We fit several galactic halo models of the stellar distribution 
to our eight lines of sight, and explore the structural parameters 
resulting from the best fits, as well as the influence of substruc¬ 
ture in those parameters. We find that the halo is best represented 
by a broken power law with index n„, = -2.50 ± 0.04 in the in¬ 
ner halo (R < Rbreak = 19.5 ± 0.04) and n out = -4.85 ± 0.04 
in the outer halo. Our data cannot constrain whether a change 
in the polar axis ratio also accompanies the break in the halo. 
The best fit values for the polar axes ratio indicate a quite 
oblate halo: q = 0.63 + 0.02. The simpler (non-broken) triaxial 
power law models favour a practically axisymmetric halo, with 
w > 0.81 + 0.07 and the rest of parameters equal to those of the 
axisymmetric one. 


We find that fitting models to data that contains substantial 
substructure can bias significantly the perception of triaxiality, 
decreasing the disk axis ratio w by 15%. We also find that dif¬ 
ferent distance modulus bin sizes and the inclusion or exclusion 
of particular lines of sight can moderately influence our mea¬ 
surements of some structural parameters. This calls for carefully 
crafted analysis and tailored tests in any future studies. When 
compared to previous works, the choice of stellar tracer seems 
to have no significant influence on the values of the structural 
parameters, at least for these distance ranges. 


Comparing our density profiles to the smooth model fits, we 
recover the presence of the Sagittarius stream in groups G and 
H. The Sagittarius stream in the direction of group H seems to 
extend further out from the Galactic centre than the models have 
so far predicted, and confirms previous RRLyrae detections as- 
sociated with the stre am at such distances ( Ibataet akl (1200 ll) . 
iTotten & Irwinl (119981) and lDohm-Palmer et alJ (1 2001 If ). We also 
find evidence of more modest substructures extending over a 


long range of distances in group C (35 < Rgc ^ 60 kpc) and 
quite concentrated in distance in group D (Rcc ~ 20 kpc). 

Our pencil beam survey has demonstrated that even a rela¬ 
tively small numbers of narrow fields of view, provided they are 
sampled sufficiently deep and with an abundant tracer, can place 
competitive limits on the global density profile and shape of the 
Galactic halo. The advent of similarly deep, wide-area surveys 
-like KiDS, VIKING and LSST- therefore promises to enhance 
substantially our understanding of the halo. 
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